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ABSTRACT 

In a previous paper [1], a new merit function for determining the strength 
performance of flawed composite laminates was presented. This previous analysis 
was restricted to circular hole flaws that were large enough that failure could be 
predicted using the laminate stress concentration factor. In this paper, the merit 
function is expanded to include the flaw cases of an arbitrary size circular hole or 
a center crack. Failure prediction for these cases is determined using the point 
stress criterion. An example application of the merit function is included for a 
wide range of graphite/epoxy laminates. 


INTRODUCTION 

In a previous paper [1], a new merit function was determined to examine the flaw 
tolerance of composite laminates. In brief summary, it was determined that the 
strength performance of a uniaxially loaded infinite width orthotropic sheet could 
be optimized by maximizing the new merit function, yfE x . In this function. Ex 
is the longitudinal extension modulus of the laminate and \|/ is the flaw tolerance 
factor used to predict the reduction in strength of a flawed laminate. The values 
of \| / are < 1 and can be used to account for various flaw effects such as holes, 
cracks, or impact damage. For circular hole flaws that were large enough that 
failure could be predicted using the laminate stress concentration factor, it was 
shown in [1] that if/ was equal to 1/k, where k is the orthotropic stress 
concentration factor. Therefore, the merit function for the case of a large 
circular hole is Ex/k. For a small circular hole or a center crack, since failure 
cannot be predicted through the use of a simple stress concentration factor, the 
same merit function cannot be used. 

In this paper, the merit function is expanded to include the flaw cases of an 
arbitrary size hole or a center crack. In order to predict the failure strength of 
an infinite width plate in the presence of these flaws, the point stress criterion [3] 
is used. After development of the merit function for each of the flaw cases, 


1 



examples are presented of the application of the merit function to a wide range of 
graphite/epoxy laminates. 

DESIGN OF A LAMINATE FOR IMPROVED FLAW TOLERANCE 
IN THE PRESENCE OF A CIRCULAR HOLE 


For an infinite width orthotropic sheet containing a hole of radius R, if a uniaxial 
stress a is applied parallel to the y-axis (Figure 1), then the normal stress, dy, 
along the x-axis is approximated [5] by 


(J y (X,0) = 




+ 3 




( 1 ) 


where k is the orthotropic stress concentration factor for an infinite 
orthotropic plate [2] 


k — 1 + 




width 


(2) 


The subscripts 1 and 2 denote the material principle axis. 

The stress distribution given by (1) is actually a modification of the isotropic 
stress distribution solution. For isotropic materials, the stress concentration 
factor, k, equals 3 and expression (1) is exact. Since (1) is an approximation of 
the stress distribution in front of a hole for anisotropic materials, the accuracy of 
this expression was tested before use. The approximate stress distribution was 
compared with the exact stress distribution solution from Savin [4]. The details 
of this comparison are contained in Appendix A. 

In order to predict failure in the presence of a circular hole, the point stress 
failure criterion [3] has been used. Using the point stress criterion, failure of a 
flawed material is assumed to occur when the normal stress, G y, at a certain 
distance, do, ahead of the flaw reaches the ultimate strength of the unflawed 
laminate, Go, or when 


< 3 > 

do is a material parameter that needs to be experimentally determined for each 
laminate under consideration. 
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Using the point stress criterion along with the stress distribution given by (1) 
results in the flawed to unflawed strength ratio for an infinite width plate 
containing a circular hole as 


where 


= 2 

{2 + ^ 2 +3^ 4 -(/c-3)(5| i 6 -7| 1 8 )} 



R 

{R + d 0 ) 


(4) 

(5) 


In (4), <r^ is the ultimate strength of the flawed laminate. 

Using the methodology presented in [1], expression (4) is used to define the flaw 
tolerance factor for the flaw case of a circular hole, y / CH , as 

VCH ” {2 + £, 2 +3^ -(k - 3)(5tf - 7^, 8 )} ' (6) 

As R — > 0, then i/r cw — > 1 which, as expected, implies that there will be no 
reduction in strength. As R becomes large ( > 2.5 in), then \f/ CH — > 1/k. This is 
shown in Figure 2 where the value y/ CH *k is plotted against hole radius for three 
values of the stress concentration factor. Since \f/ CH — > 1/k, then \f/ CH *k — » 1 . 
For a hole diameter larger than 5 in, each of the curves shown in Figure 2 
approach this limit with the error only reaching a maximum of 1 1 % for the k=5 
curve. Since expression (6) returns the expected limits of the strength reduction, 
it is used herein for a circular hole of any size. 

For the flaw case of a circular hole, the merit function is 

2E 

¥CH E * ” f + £ 2 + 3£ 4 - - 3)(5£f - 7^| 8 )| ’ <7> 


DESIGN OF A LAMINATE FOR IMPROVED FLAW TOLERANCE 
IN THE PRESENCE OF A CENTER CRACK 

For a center crack of length 2c in an infinite width orthotropic sheet under an 
axial load a (Figure 1), the exact anisotropic solution for the normal stress, Uy, 
in front of the crack [2] is 


3 



ax 


( 8 ) 


a y (x,0) = 



-c 


2 


Using the point stress criterion along with (8) results in the flawed to unflawed 
strength ratio for the flaw case of a center crack [3] as 


where 



(9) 

( 10 ) 


Using (9), the flaw tolerance factor for the flaw case of a center crack, y/ C c-> can 
be defined as 

'l>cc=^|i :: t!■ ( 11 ) 

It should be noted from (11) that \f/ cc is independent of laminate properties and 
is only a function of the crack size and do- This implies that for a given crack 
size and do, Ycc ls a constant- Therefore, in the case of a center crack, the merit 
function will be an optimum by maximizing the laminate property E x . 


The merit function for the flaw case of a center crack is 

¥ccE x =(^!)e x . ( 12 ) 


Example Application of Merit Function 

For the example application of the merit function, the same baseline laminate 
identified in Table 1 of reference [1] was used. To show the application of the 
merit function, the baseline laminate was perturbed from the initial configuration 
of [±45/02/±45/02/±45/0/90]2s to a different configuration. The resulting 
laminate properties and flaw tolerance factor were then determined for each of 
the resulting laminates. 

To do this, the baseline laminate was considered to consist of only three layers. 
These three layers each have thickness t t and consist of all the 0° plies, all the 
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I±45°l plies, and all the 90° plies respectively. Using t t =t t / h, h = laminate 
thickness, the relative values of the f ; 's were changed, and the angle of the ±45° 
ply was changed to some angle ±9. The laminate stiffness properties were then 
computed using equations (3.1) through (3.5) from [1] and the stiffness constants 
defined in Appendix A of [1] also. The resulting flaw tolerance and merit 
function for circular hole flaws was then determined using equations (2), (5), (6) 
and (7). For center crack flaws, equations (10) and (12) were used. 

For each of the flaw types, the dimensionless thickness f 90 was left constant at 
t 9Q = .083, its value in the baseline laminate. The values of t Q and t e were then 
varied from 0 to .917 while holding the sum t 0 +t e =.9\l a constant, and 9 was 
varied from 20° to 60°. The result these changes have on the laminate stiffness 
E x is shown in Figure 3. As expected, E x increases with increasing t 0 and 
decreasing 9. The resulting changes in the individual merit functions are 
discussed in the following sections. A value of do = .04 inches was used in all 
calculations since this value has been shown to give good results for 
graphite/epoxy laminates [3,6]. 

Circular Holes 

Three hole sizes were used: R= .125 in., R=.5 in. and R= 2.5 in.. The R=2.5 in. 
hole size is the same as that presented in [1]. The flaw tolerance factors for 
0=±45 and variable t 0 are shown in Figure 4, while the effect of varying 9 is 
shown in Figures 5, 6 and 7. A major difference between the small and large 
hole cases is the fact that the flaw tolerance factor increases with increasing t 0 for 
the small hole case (R=.125), while it decreases with increasing t 0 for the larger 
hole cases. In addition, decreasing 9 improves the flaw tolerance factor for 
small holes, while worsening the flaw tolerance factor for larger holes. The 
resulting values of the merit function are shown in Figures 8, 9 and 10. 

For the R=.125 case (Figure 8), it can be seen that the optimum values of the 
merit function will occur when f 0 is increased to its maximum value of .917 and 
when 0 is decreased to 20°. Significant increases in the merit function compared 
to the baseline laminate are possible. 

For the R=.5 case (Figure 9), optimum values for the merit function also occur at 
high values of f 0 and low values of 0, but not at the extremes of these 
parameters. Significant increases in the merit function over the baseline laminate 
are still possible. 
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The results for the R=2.5 case (Figure 10) are identical to those presented in [1], 
Optimum values of the merit function now occur at various combinations of t 0 
and 0. Increases in the merit function are still possible, but are not as significant 
as in the smaller hole cases. 

Center Cracks 

Three different crack sizes were examined: cracks with half lengths of c=.125 
in., c=.5 in and c= 2.5 in.. Since the values of \\f are a constant for any given 
crack size, the optimum values of the merit function will occur when E x is 
maximized. The values of \j/ for each case are shown in Figure 4. The resulting 
merit function values are shown in Figures 11, 12 and 13. As intuitively 
expected, larger cracks reduce the strength more than smaller cracks. The merit 
function for each crack size is optimum at high values of / 0 and lower values of 
0 since these values maximize E x . In each case, significant increases in the merit 
function compared to the baseline laminate are possible. 


Discussion 

The example application of the merit function uses the same graphite-epoxy 
laminate as presented in reference [ 1 ] . In addition to large holes, small holes and 
center cracks were studied in this example. From this example it was seen that: 

1) For the flaw case of a small hole (R< .125 in), significant increases in the 
merit function over the baseline laminate are possible. The large increases are 
driven by the fact that increasing the longitudinal stiffness E x also increases the 
flaw tolerance factor for the small hole case. For the larger hole cases, V CH 
decreases with increasing E x and the increases in the merit function are not as 
significant. 

2) For the case of a center crack, the merit function is optimized by increasing 
the longitudinal stiffness E x of the laminate since the flaw tolerance factor is a 
constant not dependent on laminate properties. Large cracks reduce the laminate 
strength more than small cracks. 

Increases in the merit function are again accompanied by significant decreases in 
the shear stiffness G X y and fluctuations in the major Poisson’s ratio as detailed in 
[1]. These changes may not be acceptable depending upon the application. 
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As detailed in Appendix A, values of the merit function at high values of / 0 may 
not be accurate for the circular hole flaws. This is due to the fact that the 
accuracy of the approximate stress distribution (equation 1) is not very good for 
highly fiber-dominated laminates. Since the merit function derived for center 
cracks uses the exact stress distribution (equation 8), the values of the merit 
function at high t 0 values should still be accurate in the center crack examples. 

It should also be noted that the results presented herein are for infinite width 
plates only. The effects of finite width are beyond the scope of this paper, and 
should be examined further. 


CONCLUDING REMARKS 

In this paper, the new merit function \|/E X developed in [1] was expanded to 
include the flaw cases of arbitrary size circular holes and center cracks in an 
infinite width orthotropic sheet. In reference [1], the merit function was limited 
to holes large enough that failure could be predicted using the laminate stress 
concentration factor. Since failure of a plate containing a small hole or a center 
crack cannot be predicted through the use of a simple stress concentration factor, 
a different merit function must be used. 

In the present paper, the merit function for the flaw cases of an arbitrary size 
circular hole or a center crack were derived using the point stress failure 
criterion. An example of the use of the merit function was presented for each of 
the flaw cases using a wide range of graphite/epoxy laminates. 

From this example, it was shown that significant increases in the merit function 
compared to a baseline laminate are possible. The increases are tempered by the 
flaw size, with larger flaw sizes having smaller increases in the merit function. 
For center crack flaws, increases in the merit function are governed by increases 
in the laminate longitudinal stiffness E x only and not the flaw tolerance factor. 
This is due to the fact that the flaw tolerance factor for center cracks is a constant 
not dependent on laminate properties. Accompanying the increases in the merit 
function are corresponding decreases in the laminate shear stiffness and changes 
in the major Poisson's ratio. These effects are fully detailed in reference [1], 
The effect of finite width on the results contained herein needs to be further 
examined. 
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ictor approaches 1/k with increasing hole radius. 
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Figure 3 - Extensional Modulus for example laminate. 




Figure 4 - Behaviour of the flaw tolerance factor for circular hole flaws and center cracks. 
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Figure 5 - Flaw tolerance factor for a flaw width of .25”. 
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Figure 7 - Flaw tolerance factor for a flaw width of 5.0”. 
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Figure 8 - Merit function for example laminate, hole radius =.125 




Figure 9 - Merit function for example laminate, hole radius 
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Figure 10 - Merit function for example laminate, hole radius = 2.5”. 




Figure 1 1 - Merit function for example laminate, crack half- width = .125”. 
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Figure 1 2 - Merit function for example laminate, crack half-width 
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Figure 13 - Merit function for example laminate, crack half-width = 2.5 




Appendix A: Comparison of the Approximate and Exact Stress 

Distribution Solutions for an Orthotropic Plate. 

In order to test the accuracy of the approximate stress distribution solution shown 
in equation (1) of this paper, the solution obtained using the approximate 
expression was compared with the exact stress distribution obtained using the 
complex variable mapping approach by G. Savin [4], The normal stress 
distribution in front of a circular hole in an infinite orthotropic plate under an 
axial load o was computed for a broad range of laminates using both methods, 
and the error between the solutions was determined. 

The exact normal stress distribution, < 7 y , along the x-axis in front of a circular 

cutout of radius R in an infinite width anisotropic plate under an axial load <7 
(Figure 1) is given by 


°»K-°) , , n c 

1 

-S 2 (l-iS 2 ) 

_| 

S,(l-iS,) 1 

a 

s,s 2 

V^ 2 -|-Sf 


1 

1 i 

k 2 -l-Sl 

^ + ~ 1 -S 2) 


(1A) 


where 


( 2A > 

K 

For an infinite width orthotropic plate. Si and S 2 are the roots of 

s4 + (g’~ 2v 0 s2 + § L=0 (3A) 

where the subscripts 1 and 2 denote the material principle axes which are aligned 
along x and y respectively. The roots of (3A) can fall into the following three 

cases: 

1) S 1 =i6,,S 2 =i62,Bi62>0. 

2) Si=S 2 =iB, B>0. 

3) S|=a+ iB, S 2 =-oc+ iB, a,B>0. 

These roots are for the case of the principle axis parallel to the x-axis and the 
loading normal to the x-axis. To rotate the principle axis by 90°, the original 
roots S| and S 2 should be replaced by 1/S 1 and 1/S 2 . For the three root cases the 
following substitutions need to be made: 
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Case 1) replace Bi and B 2 by 1/Bj and I/B 2 . 

Case 2) replace B by 1/B. 

Case 3) replace aby a/(a 2 +B 2 ), B by B/(a 2 +B 2 ). 

In order to compute the exact stress distribution, the software package 
Mathematica was used. After first finding the roots of (3 A) for the 16 different 
laminates shown in Table 1A, the exact and approximate stress distributions were 
computed and compared for each of the laminates. The material properties for 
the laminates were obtained from references [6] and [7], The stress distributions 
calculated using both methods and the error between the two solutions are shown 
in Figures 1A through 16A. In all of the figures, the dashed line represents the 
approximate stress distribution and the solid line is the exact stress distribution. 

The results of this investigation show that the approximate expression does 
provide very accurate solutions except in a few laminate cases. These cases 
include highly fiber-dominated laminates ( mostly 0° plies) or when the laminate 
consists mostly of ±45° plies. Yet even in ±45 laminates (i.e. Poe 11, 15) the 
error between the exact and approximate solutions is quite low (~6%). Since 
these laminates are the extremes of laminate design, their use is fairly limited. 
For the broader general class of laminates, the approximate stress distribution is 
accurate enough and has been used in this paper. 

The desire behind using the approximate expression given by equation (1) is that 
it is a simple closed form solution for the stress distribution that can be easily 
used without having to resort to specialized mathematical software such as 
Mathematica. While obtaining the exact solutions has become much easier with 
the aid of such packages, the exact solution technique is still very cumbersome 
and the use of simpler formulas greatly simplifies the analysis. 
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Laminate 

Ev(GPa) 

E v (GPa) 

G xv (GPa) 

EHW 

Rhodes Orthotropic (Ref. 6) 
[±45/0 2 /±45/02/±45/0/90]2s 

70.2 

34.8 

20.2 

.48 

Rhodes Quasi-Isotropic | 

(Ref. 6) 

53.3 

53.3 

20.2 

.32 

Poe 1 (Ref. 7) 
[0]8T 

129 

10.9 

5.65 

.312 

Poe 2 
[02/90/0 ] s 

100 

40.7 

5.65 

.0836 

Poe 3 

[0 2 /45/02/-45/0 2 ]s 

103 

17.9 

12.7 

.551 

Poe 4 

[90/0] 2s 

70.5 

70.5 

5.65 

.0482 

Poe 5 

[90/0/90/0/45/0/-45/0]s 

76.4 

47.4 

12.7 

.214 

Poe 6 

[45/0/-45/0] s 

75.3 

23.3 

19.7 

.649 

Poe 8 

[±45/0/±45/0] s 

50.0 

25.6 

26.0 

.698 

Poe 9 

[45/0/-45/90] s 

51.4 

51.4 

19.7 

.307 

Poe 10 

[90/45/90/-45] s 

23.3 

75.3 

19.7 

.201 

Poe 11 

[±45] 2s 

19.6 

19.6 

33.7 

.735 

Poe 12 
[0]8T 

109 

8.32 

4.82 

.314 

Poe 13 

[45/0/-45/0] 2s 

63.4 

19.2 

16.5 

.654 

Poe 14 

[45/0/-45/90] 2s 

43.1 

43.1 

16.5 

.303 

Poe 15 

[±45] 2s 

16.7 

16.7 

28.2 

.730 


Table 1A: Laminates used to investigate the accuracy of the 

approximate solution. 
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Exact and Approximate Stress Distributions 


[+- 45 / 02 / + - 45 / 02 / +- 45 / 0 / 90 ] 2 . 
Exact 


Approximat 




Exact and Approximate Stress Distributions 


[+- 60 / 0 /-+ 60 / 0 ] as 
Exact 


Approximate' 






Exact and Approximate Stress Distributions 


[02/90/0] s 
Exact — 


Approximate 






Exact and Approximate Stress Distributions 


[90/0] 2s 
Exact 


Approximate 




Exact and Approximate Stress Distributions 


[90/0/90/0/ 45/0/-45/0] s 
Exact 


Approximate 





Figure 8A: Comparison of Exact and Approximate Solutions for Poe Laminate 6. 





Exact and Approximate Stress Distributions 


[H — 45/0/ +- 45/0] s 
Exact 


Approximate 





Figure 10A: Comparison of Exact and Approximate Solutions for Poe Laminate 9. 





Figure 11A: Comparison of Exact and Approximate Solutions for Poe Laminate 10. 









Figure 13A: Comparison of Exact and Approximate Solutions for Poe Laminate 12. 






Figure 14A: Comparison of Exact and Approximate Solutions for Poe Laminate 13. 





Figure 15A: Comparison of Exact and Approximate Solutions for Poe Laminate 14. 






Figure 16A: Comparison of Exact and Approximate Solutions for Poe Laminate 15. 







